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CHAPTER I - INTRODUCTION 


I.1 BACKGROUND 

During the several years spent in designing a submarine, a 
great deal of effort is expended trying to improve the design so that 
the best possible vessel is eventually constructed. During the de- 
sign phases the designer will frequently examine problems that were 
found in older submarines and try to determine the cause of each 
problem so that they can be eliminated in the new design. At several 
stages during the design, computer models are used to assist in such 
things as structural design, equilibrium calculations, and internal 
arrangement. As the design progresses, physical models of the 
proposed vessel are built and tested in a towing tank to measure the 
hydrodynamic coefficients. Finally, when it is possible to deal with 
the vessel as a whole, a dynamic analysis is conducted to gauge how 
the submarine will perform in an underwater environment when all six 
degrees of freedom are available. 

While it is possible to gather dynamic information, such as the 
hydrodynamic coefficients, from physical model tests, the models are 


necessarily too restricted in their motion to permit a complete 
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analysis. It iS necessary to construct a mathematical model of the 


submarine to simulate the submarine's motion and thereby obtain 
enough data for a complete analysis. The simulation model will use 
the hydrodynamic coefficients obtained from the phy$ical model as a 
basis for the simulation. A properly constructed model will permit 
the designer to simulate any conceivable maneuver and gauge the 
submarine's response. The simulation mode] then becomes a tool to 
assist in improving the design and achieving the best possible 
vessel. After the submarine is designed, the model can continue to 
be useful by assisting in the evaluation of operating procedures. 
The model can be placed in any maneuvering situation without hazard 
to crew or vessel. This is especially useful in evaluating casualty 
Situations. The model can also be used to estimate the effect of 
proposed design changes to the submarine. For instance, the model 
can show the effect of changing the maximum deflection of a control 
surface or its rate of operation. 

All of the few six degree of freedom submarine simulation models 
in existence are too expensive for daily use in the design office. 
Generally, the models are complicated to use and sometimes difficult 
to obtain. These problems nave created the need for a program that 
is both inexpensive enough to permit daily use in the design office 
and simple enough so that anyone can use it. This need has formed 


the motivation for this thesis. 


1.2 MODEL DEVELOPMENT 


A mathematical model must use Newton's Law of Motion and the 





differential equations resulting from a dynamic analysis of the 
submarine. А submarine has movable appendages, so these must also 
be accounted for in the model.  Newton's Law of Motion can be 


expressed as: 


(1) Force = d (momentum) 
(2) Moment - dt (angular momentum ) 


If a submarine with a coordinate system such as shown in Figure ] is 
used, then equation (1) can be applied along each of the three axes. 
Similarly, equation (2) can be applied around each axis. This gives 
a total of six equations to represent the six degrees of freedom of 
the submarine. These six equations of motion are well-known [1] so 
they are not developed in detail here. They are, however, included 
for reference in Appendix А. | 

While the equations of motion for the submarine form the heart 
of the model, they are insufficient by themselves. On an actual 
submarine, the officer of the deck orders the rudder or dive planes 
moved in order to maneuver the ship. The routines will normally 
sense where the ship currently is and where it is supposed to be and 
then move the appendages in the appropriate direction just as the 
deck officer would do. The six equations of motion together with 
the appendage control subroutines constitute the vital components of 
the simulation model. А main program is necessary to coordinate the 
input, output, and to control the action of the subroutines. 

The components of the model are discussed in detail in Chapter 
II of this thesis. The model is written as a computer program whose 


features are discussed in Chapter III. The final chapter will 
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program will be provided in Appendix B. 


A listing of the 
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CHAPTER II - THE SIMULATION MODEL 


II.1 GENERAL CAPABILITIES 

The simulation model developed in this thesis will provide 
trajectory information for routine submarine maneuvers such as 
turning or changing depth. The trajectory, which is referenced to 
a fixed coordinate system,.is computed as a function of time. As 
the trajectory of the model is developed, the velocities and 
accelerations are calculated and stored for output to the user. 

Both angular and linear velocities and accelerations are provided. 
The angle of deflection of each control surface appendage is computed 
for each step of the entire trajectory or maneuver. 

The rate of deflection of each control surface and its angle 
of maximum deflection are specified by the user. Any of the control 
surfaces can be "jammed" by specifying a particular angle of 
deflection. The model can be initially placed on any course and 
depth and then ordered to come to any new course and depth. The 
user can select the initial speed and the appropriate thrust 


coefficients to cause a change in speed. 
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II.2 COORDINATE SYSTEMS 


The coordinate system plan used in the simulation model is 
based on [1]. The model uses two orthogonal coordinate systems; one 
remaining fixed at the water surface while the othef travels with the 
Submarine to act as a local reference system. The fixed system, 
designated Хо» Уб» 20: is related to the moving system, designated 
X, y, z, through the angles Ф, 9, v as illustrated in Figure 2. 
Quantities measured in the moving system can be referenced to the 
fixed system by using the transformation [2] given in Appendix A. 

The origin of the moving system is at the center of gravity of 
the submarine. This has the advantage that only the principal 


І зІ 0). 


yz xr 
It has the additional advantage that the equations of motion in [1] 


moments of inertia, - Lys І.» are non-zero (i.e.: 
also use this reference point for the x, y, z system. It would have 
been possible to use the centerline of the submarine for the origin. 
This would have the advantage of making better use of the symmetry 
of the vessel, however it would have the disadvantage of having to 
correct both the equations of motion and some of the hydrodynamic 
coefficients for the new origin. 

Appendage movements are measured with respect to the moving 
coordinate system. All velocities and accelerations are measured 


along the axes of the moving system. This is also shown in Figure 2. 


I1.3 EQUATIONS OF MOTION 


The general nature of the six degree.of freedom model requires 


the use of the six equations in their non-linear form. A large 








14 


c eJnbtJj 


| suuəuuolm * ‘s2204 'sen20l9A 'salbuy "sexy jo diusuonelss 
А 


PA 
18” en аа, 


схмлүан 09113 


х 


d‘y 





J'N р Эг | 
қат т № - | | 
` = Е 
\/ Ms 





15 
number of hydrodynamic coefficients are necessary for this model. 


It would be desirable to derive each coefficient from existing hydro- 
dynamic theory. Some coefficients have been obtained for bodies of 
revolution or other mathematically amenable shapes; however, when 
appendages such as control surfaces, fairwaters, and propellers are 
taken into account, the accuracy of the theoretical values is brought 
into question. For this reason it is standard practice to obtain the 
numerical value of the coefficients by experimental means. This 
usually entails the towing and measurement of physical models. 

Once the hydrodynamic coefficients are known, the equations 
can be solved for the accelerations. The six equations of motion 
must be written such that the highest order derivatives, namely 
u, v. W, D, q, r and their coefficients, appear on the left hand 
side of the equation. This gives each equation a form like: 


u + y ta. ү + а. 


Ны іры i,j*2 i,je3 P * 94,544 9 7 AG jes " 
= ғ. (u, у, м, р, Ф, г, Ф, 0, У» бү? бу? бе» I My e e ) 


The left hand side is placed in a matrix format. The six equations 


are then represented as: 


Ш 


“242902 »< . ° 
=h ocho cho oho 7h «ch 
о (бп + Ç N ~ 
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The method of solution is an iterative technique in which initial 


values are used in the right hand side functions and the six acceler- 
ations are solved for on the left side. The accelerations are then 
used to update the right side functions. A small time increment is 
made and the accelerations are again solved by using the latest 
value for the right side functions. 

With each iteration the accelerations are used in a Taylor 
series expansion to calculate tne velocities. The calculations 


have the form: 


u u 7 
v y v 
"N (tat) d" N (t) + (да) 
p с р 
4 4 4 
r r n 


The pitch, roll, and yaw angles are calculated in a similar manner 


using the angular velocity and acceleration. 


$ $ P ) At 
0 (t +At) = 0 (t) + чаъ° (4%) + q co 


V y r r 

In order to plot the trajectory of the submarine, a coordinate 
system transformation is performed to obtain the linear velocities 
23 Vos 2 in the fixed coordinate system. The velocities are used 
in a Taylor series expansion to obtain the position of the center of 


gravity of the submarine in the Хо» Ус» Žo system. 
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О О О 
у (еа = у (0) + Лу ` (at) 
% % % 


When the position of the submarine is calculated the time is advanced 
one increment and another iteration begins. 

The iteration loop will continue to operate until some test 
criteria are met. For instance, if the model is performing a dive, 
then the model tests the value of 22 (0 see if the model is at the 
appropriate depth. А test is also made of the pitch angle to see if 
it is within some specified range of zero. Lastly a check is made to 


ensure that all of the dive planes are at zero deflection. 


11.4 CONTROL SURFACES 

The motion of the model is controlled by the movement of the 
control surfaces. The control surfaces include the rudder, the 
stern planes, and the sail planes. Рог analysis of submarine designs 
the movement of the control surfaces can be governed by an automatic 
control system. In the model each control surface is provided with 
its own automatic control system. The principle of operation for 
all of the automatic control systems is the same. In each case the 
deflection of the control surface is made proportional to an error 
signal and the rate of change of the error. The sail planes, for 
instance, control the depth of the model. The calculated deflection 
of the planes is proportional to the depth error and the rate of 


change of depth. 
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d E К, (present depth - ordered depth) 


C 

+ K (rate of change of depth) 
A large error will initially cause a large deflection, but as the 
rate of change increases, the deflection will decreáse until some 
moderate rate is achieved. Since movement of the sailplanes does 
not have any significant effect on any motion other than depth, 
its control system uses only the variables associated with depth. 

The stern planes control pitch angle as well as depth. The 
control system accounts for this by using the pitch angle, g, and the 
rate of change of 9. 

5 = К: (present depth - ordered depth) 

+ Kg (rate of change of depth) 
+ K; (pitch angle) 
+ К (rate of change of pitch angle) 

The K's in the control systems are proportionality constants 
known as gain. The value of the gain determines the sensitivity and 
response of the control system. Since the error signal continually 
varies during a maneuver, it is best to keep the gain low. 1% 15 
desirable to have just enough gain on the error signal to get the 
control system moving and keep it moving in the right direction. 


The gain on tne rate signals should be much stronger to dampen out 


the oscillations caused by response to the error signal. 


The control surfaces move at a rate specified by the 
user. Whenever the calculated deflection differs from the 
present deflection, the control surface will move toward 


the calculated value at the specified rate. 
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CHAPTER THREE - PROGRAM USER'S GUIDE 


111.1 GENERAL FEATURES 

The simulation model consists of a main program and four 
subroutines. The four subroutines are named FUNC, RUDDER, DEPTH, 
and STERN. The main program reads the input data, changes units, 
initializes values, and prints the output. The main program performs 
the Taylor series expansions, the coordinate system transforms, and 
provides the logical statements for calling the subroutines. On the 
first iteration, when time equals zero, no subroutines are called; 
the output from this first iteration is then a statement of the 
initial conditions of the problem. Each succeeding iteration will 
call subroutine FUNC to calculate the new position of the model. 
Calls to the control surface subroutines are made on the basis of 
need, with no subroutine being called more than once in the iteration. 
The decision on whether or not to call RUDDER, DEPTH, or STERN is 


contained in a series of logical IF statements in the main program. 


111.2 SUBROUTINE FUNC 


Subroutine FUNC contains the six equations of motion. The 
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input parameters to FUNC consist of the hydrodynamic coefficients, 


the model velocities, the orientation in space, all of the ship's 
characteristics such as length and moments of inertia, the deflection 
of each control surface, and the propulsive coefficients. The sub- 
routine returns the value of the six accelerations: UDT, VDT, WDT, 
PDT, QDT, RDT. The input and output for FUNC are all contained in 
COMMON /FOUR/ and COMMON /FIVE/. The solution to the matrix equation 
in subroutine FUNC is made possible by a call to the library function 
LEQTIF [3]. LEQTIF performs a Gaussian reduction for the matrix 
equation. Any similar Gaussian reduction could be substituted if 
LEQTIF is not available. An explanation of the parameters used in 


the call to LEQTIF is found in Appendix B. 


III.3 SUBROUTINE RUDDER 
Subroutine RUDDER controls all horizontal motion for the 

model. The input parameters include: 89, К10, 15, 16, 111, 112, 
T17, T18, TLAGR, RRATE, RUDAMT, COURSE, DELR, DELT, R, PSI. The sub- 
routine returns a new value for DELR, the rudder deflection, based on 
a calculation using its present and ordered headings. The inputs are 
all contained in COMMON /THREE/ and COMMON /FIVE/. Subroutine RUDDER 
is called whenever the model is not on the desired course or when the 


rudder is deflected. 


III.4 SUBROUTINE DEPTH 
subroutine DEPTH controls the forward set of dive planes. 1% 


can control either bow planes or sail planes depending on what the 
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submarine is fitted with. This subroutine is sensitive to the depth 


error and the rate of change of depth. Subroutine DEPTH will move 
the forward planes in the direction necessary to bring the depth 
error to zero. As input parameters, the subroutine uses: Kl, K2 
Tl, T2, T13, T14, T15, T16, TEAGB, DELT, DIFF, ADIFF, ZDT, ATHETA, 
MAXANG, DCRIT. The subroutine returns a new value for DELB, the 
bow plane deflection, after each call. 

If MAXANG, the maximum dive/ascent angle, has been exceeded, 
then DEPTH will return a diagnostic write statement. The user may 
specity MAXANG, the maximum pitch angle for the model. Whenever 
MAXANG is exceeded the diving planes are moved to reduce the pitch 
angle and a diagnostic signal is generated from within subroutine 
DEPTH. The diagnostic will say “Maximum dive/ascent angle exceeded 
at time" .. Standard fixup taken." Since the dive planes only 
begin to react when MAXANG is exceeded, the submarine will overshoot 
the angle before a reduction in the angle actually occurs. The 
diving planes will continue to operate until the pitch angle is 
less than MAXANG. 

DEPTH makes only a minimal attempt to control the pitch angle 
of the model; it only warns the user when the angle is exceeded, and 
then moves the bow planes so that the pitch angle does not grow 
larger. The principal purpose of DEPTH is to provide depth control; 


it does this in conjunction with subroutine STERN. 


111.5 SUBROUTINE STERN 


Subroutine STERN controls the movement of the stern planes to 
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achieve the desired depth and pitch angle. The movement of the 


planes is sensitive to the depth error, the rate of change of depth, 
the pitch angle, and the rate of change of pitch. The stern planes 
will move to bring the depth error to zero and the pitch angle to 
zero. As input parameters the subroutine uses: 65, K6, K7, K8, T3, 
T7, T8, TIO, TLAGS, STERAT, STERMX, THETA, Q, DELS, DIFF, ADIFF, 

ZDT, ATHETA, MAXANG, NOPICH, DCRIT. The subroutine returns a new 
value for the stern plane deflection, DELS. The input and output 
parameters are all contained in COMMON /ONE/, COMMON /FIVE/, and 
COMMON /SIX/. STERN is called to achieve a depth change or to correct 


the pitch angle. 


111.6 RANGE VARIABLES 

Since it is not usually possible to bring a computer simulation 
model to a precise depth or angle, it is necessary to define ranges 
around the desired depth or angle which will be acceptable to the 
user. For example, if the submarine were to make a depth change of 
900 feet, the dive may be considered complete if the model settled 
within ten feet of the desired depth. The range is specified by the 
user to enable him to achieve whatever precision is desired. Within 
this range the main program will not make a call to the control sur- 
face subroutine; since no call is made, the control surfaces cannot 
be moved. In the program, the variable name for the depth range is 
NODIFF; it is usually set at + 5 or + 10 feet. Subroutine DEPTH will 
be called whenever the model is off of its-ordered depth by more than 


NODIFF. The range about zero pitch angle is given the name NOPICH; 
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it is usually set at * 1 degree. STERN is called if the pitch angle 


is greater than NOPICH. The variable name for achieving the proper 
course is ONCRS; it is usually set at + 1 degree. RUDDER will be 
called whenever the model is off course by more than ONCRS. 

An additional range variable is used during diving or surfacing. 
When the model moves within some critical range of the desired depth, 
it is time to begin moving the dive planes to zero angle and let the 
vessel glide into the desired depth. This maneuver prevents 
oscillation of the dive planes as the error signal and the rate 
Signal both become small. This also ensures that the planes are at 
or near zero angle by the time the desired depth is reached. The 
name of this variable is DCRIT; it is usually set at 50, 75 or 100 
feet depending on the speed of the submarine. STERN will be called 


whenever the model is off its ordered depth by more than DCRIT. 


III.7 PROGRAM PARAMETERS 

The iterative nature of the program relies on a time increment 
being made after each step. The size of the time increment is 
optional; the smaller the increment the more accurate the calculations. 
In selecting the size of the time increment, one must bear in mind 
the length of time required to complete the maneuver and the storage 
capacity of the program. Due to the quantity of information that is 
calculated by the program, it is practical to print only half of the 
data at one time. The other half of the data is stored in the 
array STOW. This data is then printed when the maneuver is complete. 


STOW has the capacity to hold the data resulting from 600 iterations. 
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A time check is incorporated in the main program to enable the 


user to limit the number of iterations. Achieving a number of 
iterations equal to the variable ICNT will cause the program to stop 
the present maneuver, print out all data, and see if another maneuver 
is desired. ICNT should not be made larger than 600 so that the 
available storage space is not exceeded. 

Each control surface subroutine has a time lag scheme which 
senses the initial command to the control surface and prevents 
immediate action. A time lag occurs each time the direction of 
movement is changed. A different time lag may be specified for 
each control surface. 

The variable INDEX is used to allow the user to perform more 
‘than one maneuver with a given submarine and then shift submarines 
to perform more maneuvers. If INDEX is less than or equal to zero 
then the same submarine coefficients will be used for each set of 
initial conditions. If INDEX is greater than zero the program will 


read new coefficients as well as a new set of initial conditions. 
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CHAPTER FOUR - TEST RESULTS AND CONCLUSIONS 


IV.1 FUNDAMENTAL MOTION TEST 

The validation of the simulation model is a matter of impor- 
tance. The method used was an independent test of each component of 
the model, and then a series of tests with the components of the 
model working together. The test maneuvers were selected either 
because the correct dynamic response was known or because the maneuver 
was simple enough so that the general nature of the response could 
be predicted. 

The first portion of the model to be tested was the subroutine 
FUNC. FUNC obtains the solution to the six equations of motion. It 
was important to establish that these equations were properly 
installed in the program. А test of subroutine FUNC necessitated a 
simultaneous test of the MAIN program to read in the hydrodynamic 
coefficients, set up the "A" matrix for FUNC, perform the Taylor 
series expansions, and write out the results. The test for FUNC 
was to reduce the GM of the vessel to zero and then to assume a 
constant angle on the dive planes. If the.model was working cor- 


rectly it should traverse a perfect circle in the vertical plane. 
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It was known that the rudder should not move and the model should 


not roll or yaw. The model should assume a circular trajectory with 
a constant angular velocity and a constant vertical component of 
velocity, w. Figure 3 shows the results of the maneuver. À circular 
trajectory was quickly achieved. The angular velocity was constant, 
w was constant, and the position of the model in the fixed coordinate 
system confirmed the circular path. The model did not roll or yaw 


and the rudder did not move. 


IV.2 HORIZONTAL MOTION TEST 

With the knowledge that the equations of motion, the Taylor 
series expansions, and the coordinate system transforms are operating 
properly, the subroutine RUDDER was the next component to test. It 
was decided that a simple left turn of 40 degrees would be an 
appropriate test. A left turn was chosen because that presents the 
greatest opportunity for error. The original course would be 000 
degrees true and the new course would be 320 degrees true. Since the 
yaw angle is positive when turning right, the model must cope with a 
negative yaw angle as well as a proper method for dealing with course 
headings given in true bearing. The subroutines DEPTH and STERN were 
rendered inoperative to give the opportunity of seeing RUDDER operate 
without interference. The model would be checked to ensure that the 
proper rudder rate was used and that the rudder angle did not exceed 
the maximum angle ordered. The model would be expected to roll and 
Squat in the turn. A change in depth was expected since the dive 


planes could not act. As the model approached the new heading the 
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rudder should be put amidships and the vessel should steady out. 


In Figure 4 the relationship between the rudder and the heading 
is shown as a function of time. The model reacted as predicted 
except possibly at the end of the turn. The rudder was amidships and 
the new heading was achieved but the program did not run long enough 
to ensure that the model was steady on the course. The model did roll, 
Squat, and change depth as predicted. The test for RUDDER was 
considered to be accurate and sufficiently complete to warrant moving 


to the next test. 


IV.3 VERTICAL MOTION TEST 

Since DEPTH and STERN operate together to regulate the depth 
and pitch of the model, they were tested together. The test consisted 
of a simple dive with a depth change of 700 feet. The course was not 
changed so the rudder should not move. The model should not roll or 
уан. The dive planes should move in the proper direction and at the 
ordered rate. They should not deflect more than the ordered angle. 
The model should pitch downward and if the maximum ordered pitch 
angle is 222000 then the dive planes should move to reduce the 
pitch angle. As it approaches the desired depth, the model should 
slow its rate of descent and settle within ten feet of the depth, 
with 0 | degree of pitch angle. 

The trajectory for this test is shown in Figure 5. The model 
achieved steady state only six feet beyond the desired depth; the 
pitch angle was -.03 degrees. The dive planes were all at zero 


deflection. The model stayed within ten feet of the ordered depth 





29 


Ор 
(Ѕәәлбәр) 


ISq 
Ou 


08: 


ша 


р әмпбі- 


SWIL JO Uonounj 


JO eibuy puy uono»eljeg чәррпу 


әш 


ISg 


E 


UOI12Ə|JƏəQ 
Јәррпу 


О 


S 


OL 


Gl 


(гоомбар) 
Ччопзецза 


4әрргу 





30 


G aynbiy 


ома әіацс y JOY ошу 10 UOUN у Sy uideo 


әш! | 





® 

009 

OOG 
(1991) 

ООЁ 
эбиецэ 

00$ 
чзаәа 


DOC 


ОО! 








3l 
and within one degree of zero pitch angle for twenty seconds. This 


was done to ensure that a steady state had been achieved. 


IV.4 COMPLETE MODEL TEST 

As a final test the entire simulation model must work together. 
This test repeated the forty degree left hand turn but this time the 
dive planes were allowed to react to try to maintain the depth. The 
program was kept running until the pitch angle was within one degree 
of being zero, the model was within one degree of the proper course 
and the depth was within ten feet of the ordered depth. The model 
was initialized on course 000 degrees true at an initial speed of 
twenty knots. 

The maneuver was successfully completed. The model remained 
. Steady at 320 2] degrees true and the final depth was within one 
foot of the ordered depth. Both the pitch and roll angles were near 
one degree and were decreasing in magnitude. All of the control 
Surfaces were at zero angle of deflection. It should be noted that 
the new course had been achieved during the first sixty seconds and 
that the course was maintained by the model while the proper depth 


and attitude were being obtained. 


IV.5 USE AS A DESIGN TOOL 

There are numerous design tasks that could be used to demonstrate 
the simulation model; a simple example is sufficient for this demon- 
stration. Suppose a designer wanted to know the dynamic effects of 


increasing the rudder rate in a turn. The designer would elect to 





32 
use a simulation model. Since he would not want the dive planes to 


interfere with the analysis, he would set NODIFF, DCRIT, and NOPICH 
at large values so that the dive planes would not move. Then he 
would set RRATE, the rudder rate, at the value he désired to test 
and order the model to come to a new course. For this run, let's 
assume that the designer performed the left hand turn from 000 
degrees true to 320 degrees true. Не used a rudder rate of two 
degrees/second. He then caused the model to perform the same 
mManuever again with a new rudder rate of four degrees/second. With 
the output from these two maneuvers he could easily find the time 
required for the turn, the advance and transfer of the submarine, 
and the roll and pitch angles as a function of time. The designer 
could use the information for whatever analysis he had in mind. The 
model could be run again-at new rudder rates or the dive planes could 
be brought into play or virtually any other maneuver could be 
Simulated. This model can also do snap roll analysis [4]. 
The cost of running the simulation model is so low that it can 
be used on a daily basis if desired. For the test turns mentioned 
above, the cost was less than five dollars, which included reading 
the cards, compilation, execution, and 900 lines of output. A 
designer who used the model regularly could have the model as an on 


line dataset which would greatly reduce the cost of using the model. 


IV.6 CONCLUSIONS AND RECOMMENDATIONS 
This simulation model does accurately simulate the six degree 


of freedom motion of a submarine for moderate maneuvers. The cost of 
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operating the model is very low, which will allow frequent usage. 


The model is designed to permit a great deal of flexibility in the 
application of the model. For those designers who use the simulation 
model, it should be a useful design tool. 

If work were continued on this simulation model, it is recom- 
mended that some time be spent in improving the appendage control 
subroutines. Further application of control theory would be helpful 
with the appendage subroutines. The data input could be organized 
more efficiently to remove some of the opportunity for error. It 
could be very useful to have a plotting routine in the model to 


visually display the information. 





(3) 


(6) 
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APPENDIX A 
А.1 NOTATION i 
Syınbol Dimensionless Form Definition 
В' = в ‚ S 
B 506203 Buoyancy force, positive upward 
CB Center of buoyancy of submarine 
са Center of mass of Submarine 
I, ; | 
Г. ге = 3242 Moment of inertia of submarine about x axis 
p Е 
І - 
P Г m Moment of inertia of submarine about y axis 
‚. _ 12 T | 
D Eu 3545 Moment of inertia of submarine about z axis 
Ix : 
E E = 205 Product of inertia about xy axis 
Iyz En. 
= [yz = ET Product of inertia about yz axes 
Izx 
Lx Lx ® 5245 Product of inertia about zx axes 
1 K . 
К = —5—; Hydrodynamic moment component about x 
dp” U axis (rolling moment) 
Ка 
K, cs 2-7; Rolling moment when body angle (m, 6) and 
рі? Ч control surface angles are zero 
i К» s 
Ren u. ma Coefficient used in representing Kx as a 
f зас U function of (n-1) 
K RU ED First order cocificient used in representing 
P Р рм K as a function of p 
К; к . : | 
K. К.' = — Coefficient used in representing Kasa function 
P P іре of: $ 
K K "а Soto Second order coefficient used in representing 
= 2 e : 5 > = 
РЇР| РЇР| ~ $4 K as a function of p 
K '= —Е- Coefficient used in representing K as a function 
Pq РЯ р of the product pq 
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Coefficient used in representing Kasa 
function of the product qr 


е 


First order coefficient used in representing 
K as a function of r 


Coefficient used in representing K as a 
function of r 


First order coeíficient used in representing 
K as a function of v 


Coefficient used in representing K as a 
function of v 


Second order coefficient used in representing 
K as a function of v 


Coefficient used in representing K as a function 
of the product vq 


Coefficient used in representing K as a function 
of the product vw . 


Coefficient used in representing K as a function 
of the product wp 


Coefficient used in representing K as a function 
of the product wr 


First order coefficient used in representing 
K as a function of 5, - 


Overall length of submarine 


Mass of submarine, including water in Їгее- 
flooding spaces 


Hydrodynamic moment component about y axis 
(pitching moment) 


Pitching moment when body angles (xy, 8) and 
control surface angles are zero 


Second order cocíficient used in representing 
M as a function of p. First order coefficient is 
zero. 


First order coefficient used in representing 
M as a function of q 


First order coefficient used in representing 
Md as a function of (7-1) 


Coeíficient used in representing Masa 
function of à 
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Second order coefficient used in representing 
M as a function of q 


е 


Coefficient used in representing М 65 аѕ а 
function q 


Coefficient used in representing Masa 
function of the product rp 


Second order coefficient used in representing 
M as a tunction of r. First order coefficient 
is zero 


Coetficient used in representing Masa 
function of the product vp 


Coefficient used in representing M as a 
function of the product vr 


Second order coefficient used in representing 
M as a function of v 


First order coefficient used in representing 
M as a function of w 


First order coefficient used in répresenting 
Мо as a function of (7-1) 


— 


Coefficient used in representing M as a function 
of w 


First order coefficient used in representing M 
as a function of w; equal to zero for symmetrical 
function 


Coefficient used in representing Mq as a function 
Of w 


Second order coeíficient used in representing 
M as a function of w 


First order coefficient used in representing 
i f (n- 

М» | as a function of (7-1) 

Second order cocfficient used in representing 

M as a function of w; egual to zero for sym- 

metrical function 


First order cocííicier.t used in representing 
M as a function of 5, 


First order coefficient used in representing 
M as a function of 6, 


First order cocfficient used in representing 
Mag s á function of (7-1) 
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N N' z PEUT: Hydrodynamic moment component about z 
axis (yawing moment) 
Nx й 
Na Na = тр Yawing moment when body angles (q@, В) and 
2р control surface anyles are zero 
N 
N N ' = Tan First order coefiicient used in representing М 
P P A as a function of p 
Ne , 
N. s Coefficient used in representing Nas a function 
P P 5р4 of р 
М. 
М N '= š Coefficient used in representing N as a function 
РЧ Ра ¿pt of the product pq 
Nar 
Nor Ыы; = ipt Coefficient used in representing N as a function 
of the product qr 
N 
N, т First order :oefficient used in representing N 
30430 as a function of r 
N 
Е 'z та First order coefficient used in representing 
p? d реу Ny as a function of (n- 1) 
М; re : : : 
N. N.' = —— Coefficient used in representing N as a function 
Е г $p of r 
,_- "Nelrl а ЕТ . ) 
М. | N Ir] :———— Second order coefficient used in representin 
5 E 324 N as a function of r 
N N 1 Мігібт ЕНТ еа а по М 
Ir! ae |r|ór E "T oetiicient used in representing 1 5r asa 
2р4 function of r 
N 
N, Ne = - First order ccefficient used in representing N 
30470 as a function cf v 
N | 
m N; '= = First order coefficient used in representing N, 
П 1 40420 as a function of (n- 1) 
N: 
N, N. ' = - Coefficient uSed in representing Nasa 
v 201 function of v 
N : 
5 = 4 Coefficient used in representing М as a iunction 
3 “4 ¿pts of the product vq 
‚_ Nivir "OM Vat : 
N | ч err m Coefficient used in representing Nyasa 
[vir ут фе function of v 
N 
N Iv] БЕЛЕ z evi Second order coefficient used in representing 
Ма spt N as a function of v 
N 
N [v| Nov] 'z un First order coefficient used in representing 
ven n 224 Nil as a tunction of (n-1) 
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Cocfficient used in representing N as a function 
of the product vw d 


Coefficient used in representing N as a function 
of the product wp 


Coefficient used in representing N as a function 
of the product wr 


First order cocfficient used in representing N 


as a function of бү 


First order coefficient used in representing 
N,, as a function of (n-1) 


Angular velocity component about x axis 
relative to flu’d (roll) 


Angular acceleration component about x axis 
relative to fluid 


Angular velocity component about y axis relative 
to fluid (pitch) 


Angular acceleration ccmponent about y axis 
relative to fluid = 


Angular veiocity component about z axis 
relative to fluid (yaw) 


Angular acceleration component about z axis 
relative to fluid 


Linear velocity of crigin of body axes relative 
to fluid 


Component of U in direction of the x axis 


Time rate of chanze of u in direction of the 
x axis 


Command spezd: steary value of ahead speed 
component u for a given proreller rpm «hen 
body angles (а, 8) and control surface angles 
are zero. Sign changes with propeller reversal 


Component of U in direction of the y axis 


Time rate of change of v in direction of the 
y axis 
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Component of U in direction of the z axis 


е 


Time trate of change of w in direction of the 


z axis 


Weight, including water in free flooding spaces 


Longitudinal body axis; also the coordinate cf a 
point relative to the origin of body axes 


The x coordinate of CB 


The x coordinate of CG 


A coordinate of the dispiacement of CG relative 
to the origin of a set of fixed axes 


Hydrodynamic force component along x axis 
(longitudinal, or axial, force) 


Second order coefficient used in representing 


X as a function of q. 


is zero 


First order coefficient 


P d 


Coefficient used in representing X as a function 


of the product rp 


Secogd order coefficient used in representing 


X as a function of r. 


Zero 


First order coefficient is 


Coefficient used in representing X as a function 


of ü 


Second order coefficient used in representing 
X as a function of u in the non-propelled case. 
First orger coetticient 1s zero 


Coefficient used in rer-esenting X as a function 


of the product vr 


Second order coeííient used in representing X 


as a function of v. First order coelficieat is zero 


First order coefficient used in representing X, ç; 


as a function of (7-1) 


Coefficient used in representing X as a function 


of the product wq 
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Second order coefficient used in representing 
X as a function of w. “First order coefficient is 
zero 


First order coefficient used in representing X. 
as a function cf (m-1) 


Second order coefficient used in representing X 
as a function of Sp: First order coefficient 
is zero 


Second order coefficient used in representing 


X as a function of бү. First order coefficient is 
zer 


First order coefficient used in representing 


Ха г as a I: nction of (7-1) 


Second order coefzicient used in representing X 
as a function of ôs. First order coefficient is 
zero 


First order coefficient used in representing 


: (11- 
Ру as a function of (7-1) : 


эт 


Lateral body axis; а150 the coordinate of a 
point relative to the origin of body axes 


The y coordinate of CB 


The y coordinate of CG 


A coordinate of th? displacement of CG relative 
to the origin of a set of fixed axes 


Hydrodynamic force component along y axis 
(lateral force} 


Lateral force when body angles (a, B) and control 
surface angles are zero 


First order coefficient used in representing 
Y as a function of p 


Coefficient used in representing Y as a function 
of p 


Second order coefficient used in representing 
Y as a function of p 
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Y 
D зоон Cocificient used in representing Y as a function 
РЧ РЧ 224% of the product pq 
Y 
Y Y SE 1- Cocfficient used in representing Y as a function 
a q pt of the product qr 
x 
n Y'= = First order coeíficient used in representing Y - 
Р 36470 as a function of r 
= ar First order coefficient used in representing 
rn гт 2200 Y, as a function of (n-1) 
Ye 
Y. Y ee Coerficient used in representing Y as a function 
E ‚г зро of È 
ЕЗІН; 
Ү Y — Cocíficient used in representing Y, asa 
[т!&г Irlör эрс U function of r = 
Y 
Y Y ' = = First order coefficient used in representing 
Ы У 20 Y as a function of v 
Y 
X = — First order coefficient used in representing 
M ыг Pot U Y as a function of (n- 1) 
Y. 
Y. Үү z М Coefficient used in representing Yasa 
М iet function of Y - 
Үт 
E Y ' s 3 Coefíicient used in representing Y asa function 
vq 4 79% of the product vq 
Y 
ir Par] E T кшн used in representing Y, as a function 
Y 
Y lv! Y | | а Зем Second order coefficient used in representing 
ыг хай 304 Y as a function of v 
Y vivim 
| Ээ : Te ! : 
Nm ЭР EU First order cocfíicient used in representing 
< И as a function of (7-1) 
' Y 
14 X = == Cocfficient used in representing Y as a 
Mad 22 pt function of the product vw 
Y 
Y 122 Cocfficient used in representing Y asa 
“Р ЖЕ 2043 function of the product wp 
Y 
X ! 24 I. Coetficient used in representing Y as a 
EE “т 20 function of the product wr 
Y 
Y ! - = : First order coefficient used in representing 
бт or 50470 Y as а function of 6r 
E : 
Ү Ү ' = Sen. First order coefficient used in represerting 
627. Srn $р4* 0 Т as a function of (7-1) 
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Normal body axis; also the coordinate of a 
point ‘relative to the origin of body axes 
, 


The z coordinate of CB 


The z coordinate of CG 


A coordinate of the displacement of СС 
relative to the origin of a set of fixed axes 


Hydrodyuamic force component along z 
axis (normal force) 


Normal force when body angles (a, B) and 
control surface angles are zero 


Second order coefficient used in representing 
Z as a unction of p. First order coefficient 
is zero 


First or ler coeíficient used in representing 
Z as a functicn of q 


First order coefficient used in representing 


Za as a funct-on of (1-1) 


a 


Coefficient used in representing Z as a 
function of q 
Coefficient used in representing Z зз а 


function of q °. 


Coefficient used in representing Z as a 
function of the product rp 


Second order coefficient used in representing 


Z as a function of r. First order coefficient 
is zero 


First order ccefficient used in representing 
Z as à function of w 


First order coefficient used in representing 


Z, asa function ot (7-1) 


Coefficient used in represcnting Z as a 
function of w 


First order coefficient used in representing 
Z as a function of w; equal to zero for sym- 
metrical function 


Coefficient ued in representing Рм asa 
function of q 
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Second order coefficient used in representing 
Z asa function of w 


First order coefficient used in representing 


E as a function of (0-1) 


Second order coefficient used in representing - 
Z asa function of w; equal to zero for sym- 
metrical function 


First order coefficient used in representing Z 
as a function оѓ ёр 


First order coefficient used in representing 
Z as a function of 6, 


First order coefficient used in representing 
Z, asa function of (7-1) 
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Angle of attack 
Angle of drift 
Deflection of bowplane or sailplane 
Deflection of rudder 
Deflection of sternplane 
The ratio — 
O — -. 
e ratio — 
Angle of pitch 
Angle of yaw 


Angle of roll 


Sets of constants uscd in the representation of 
propeller thrust in the axial equation 
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A.2 EQUATIONS CF MOTION 
AXIAL FORCE 


ЭГ - Vr + wq - xg (9° + г?) +ус (ра - r) + zG (pr +а) | = 


p M 2 2 | 
+ — DO E MUTO лэг 
2 ^ L^qq 3 rr гр ‘Р 


Еа эх ] 
вэ + Х, u +x. ea wq 


+£ ШЕ: 'u? +X 'v?4+X бад | 
2 uu vv ww 


2 
2 


2 2 72% 2 (52-2 | “а 2 
Ap ЕЛ ôr ст tX 65 95 DD са 


i 


i 2 2 2 
ue РЕ +b, uu. tcu. | 
- (W - B) sin 8 


v? +X tw? +X ô 2 u? 


p hi | 
й 2 i Кут 7 Oror] r 


242 
р Х 5557 о | (7-1) 
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WATER AL FORCE 


m[ Y - wp +ur - yG (r + p) + zG (qr - P) * xz (p +1) ] = 


| 


| + t [Yur r +Y. p + Y | +Y рась Х '' r] 
E > г Ре аа 


| EA 3 S t j ! ! | 
Е 1 ыза ұлыт Реле г wr 


1 = 
Nana, 22 
ani tv lui] 


p 3 
=, A | v. ur +Y 'up +Y |r|ór u|r|ór + Y 
i | 1,3 
12 Ее u? 4 Y, uv t IE у |(у° + w*)? | | 


+2: 42 ү ту t Y "u? бт] 
к 2 vw ór 


+ (W - B) cos 8 sin $ 


NES 
+ o 
2 t = 


п ur (1-1) 


'P ,2 i 2 2,3 | 
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NORMAL FORCE 


: . mis 
m| w- uq + vp - zg (p? * a^) + x, (тр - 4) +Ус (79 %Р) | = 


р “І 227 02222 007472 ГЭ 
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w vr vp 
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ROLLING MOMENT 
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PITCHING MOMENT 
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YAWiING MOMENT 
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AXIS TRANSFORMATIONS [2] 

1) A transformation of an axes system takes a quantity described 
in one frame of reference and transforms it into another frame of 
reference such that if we measured the same quantity in the second 
frame of reference the transformed quantity and the measured quantity 
would be identical. 

2) Transforms between frames are needed in the study of the motions 
of ocean vehicles because the equations cf motion for such a vehicle 
are most easily derived in the inertial frame attached to the earth 
(x, » Уу» 2.) frame, while the forces acting on the vehicle are most 
easily evaluated in the frame attached to the vehicle (x, y, z). 
Hence, we ultimately desire to transform the equations of motion 

from the inertial frame into the non-inertial frame fixed in the 
vehicle. 


> | > 
3) If Vo is some vector measure in the x › 2. frame and V 


о” Ўо 
Some vector measured in the x, y, z frame which is only changed 
in orientation then: | 
ү =T (4,8,9) 7 where T (v, 0,9) = the transform. 
Where 
cos@ cosy 


Т (9, 9,4) = -siny cos + sing sing cosy 


Sing siny + cos¢ cosy sing 


coso siny -sing 
cos¢ cosy + sing sing siny 511Ф С050 


-Sinó cosV + СО$Ф sing sin coso COS4 
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4) 130 | апа V are the same vectors as in (3) above, then: 
T (60,9) V 
Where 
Ц, в, У) = 
cosé@ cosy гзіпу cos? + sing sing cosy 
cosg siny cos? cos¥ + sing sing siny 
-sing зіпф со58 


sind Sind + cos¢ cosy sing 
-sing cosy + cos? sing siny 


COS С050 





B.1 LIST OF VARIABLES 


VARIABLE 


Al ,A2 


AI,BI,CI 


BORATE 


BOWMAX 


COURSE 


DECRIT 


DEL] 


DEL3 
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APPENDIX B 


MEANING 


Six by six matrix containing coefficients 
of UDT, VDT, WDT, PDT, QDT, RDT. 


Six by six matrix set equal to matrix А 
before each call to FUNC. The values of 
AA are lost in the matrix reduction performed 


by LEQTIF. 


Limits used for selecting the proper 
propeller thrust. 


Set of constants representing the propeller 
thrust in the X-equation. 


Ship's buoyancy. 

Average bowplane rate. 

Maximum ordered bowplane deflection. 
New course for the ship. 


Value of depth error when dive planes are 
returned to zero deflection. 


Calculated deflection of the bowplane. 


Calculated deflection of the stern plane. 
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нг 
За” 





EE UU UU V cS I 2 даа 


DEL4 


DELB 


DELBO 


DELGM 


DELR 


DELRO 


DELS 


DELSO 


DELT 


DIFF 


ICNT 


ID 


INDEX 
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Calculated deflection of the rudder. 


Actual bowplane deflection at time t. 


Actual bowplane deflection (units changed 
for output). 


Amount the original GM is to be changed. 
Actual rudder deflection at time t. 


Actual rudder deflection (units changed 
for output). 


Actual stern plane deflection at time t. 


Actual stern plane deflection (units 
changed for output). 


Time increment used in iteration. 


Difference between present depth and 
ordered depth. 


Six by one matrix containing solution to 
right hand side of equations of motion. 


Counter used to count the number of 
iterations 


Forty-character alphanumeric heading. 
If greater than zero, read new submarine 


coefficients. [+ less than or equal to 
zero, run same submarine for new initial 








каа ет Iu 
, P l | Я | 
2r 


Л. 
Я а 


ЛОИ РА ТА 12,12 


KCOEFF 


LEQTIF 


M 


MAXANG 


MCOEFF 


NCOEFF 


NODIFF 


NOP ICH 


ODEPTH 


ONCRS 


P,Q,R 


PDT,QDT,RDT 


PDTO,QDTO,RDTO 
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conditions. 


Moments of inertia. 


K - equation coefficients: 


Ship's overall length. 


Matrix reduction subroutine from IMSLIB [3] 


Ship's mass. 


Maximum ordered dive/ascent angle. 


M - equation coefficients. 


N - equation coefficients. 


Acceptable error range around ordered depth. 


Acceptable error range around zero pitch 
angle. 


Ordered depth. 


Acceptable error range around ordered 
course. 


Angular velocity about the x, y, and z 
axes respectively. 


Angular acceleration about the x, y, and z 
axes respectively. 


Angular acceleration about the x, y, and zZ 





PHI 


PHIO 


P0,Q0, RO 


PSI 


PSIO 


RHO 


RRATE 


RUDAMT 


STERAT 


STERMX 


STOW 


и. = 116 


ТНЕТА 
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axes, respectively (with units changed for 
output). 


Angle of roll. 


Angle of roll (with units changed for 
output). 


Angular velocity about the x, y, and z 
axes respectively (with units changed for 
output). 


Angle of yaw. 


Angle of yaw (with units changed for 
output). 


Sea water density. 

Average rudder rate. 

Maximum ordered rudder deflection. 
Average stern plane rate. 

Maximum ordered stern plane deflection. 
Storage array for half of output data. 
Present time. 

Time lag signals. 


Angle of pitch. 





THETAO 


TLAGB 


TLAGR 


TLAGS 


U,V,W 


UDT,VDT,WDT 


UO 


UOO , V0 ,WO 


WT 


X,Y,Z 


ХВ, УВ, В 


XCOEFF 


XDT,YDT,ZDT 


AGS TG, ZG 


57 
Angle of pitch (with units changed for 


output). 

Time lag for bowplane control system. 
Time lag for the rudder control system. 
Time lag for stern plane control system. 


Forward, lateral, and vertical velocities 
respectively. 


Forward, lateral, and vertical accelerations 
respectively (with units changed for output). 


Initial forward velocity. 


Forward, lateral, and vertical velocities 
respectively (with units changed for output). 


Ship's weight. 


Coordinate labels of the fixed (xo; Уу» 
21 coordinate system. 


The x, y, z position of the center of 
buoyancy. 


X - equation coefficients 


Velocities in the fixed coordinate system 


axes respectively. 


along the Ху» Yo? Zo 


The x, y, z position of the center of 











> 


mE 





YCOEFF 


ZCOEFF 


gravity. 
Y - equation coefficients 


L - equation coefficients 
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B.2 COMPUTER PROGRAM LISTING 


(z2*1-1*(1)34300N) (01*5)ay3y 
(€2*1-1*(1)44300u) (01:9)0у39У 
(4A 121*(1)4433003). (01*5)ay3y 
(22*1«1*(1)433097) (01*5)0v38 
(2241-1“(1)44302А) (01*5)av38 


(91*121*(1)343300X) (666=0N3*01*5)0v3Y4 € 


€*2*2 (X3aNI)JI 
X3GNI (666=0N3°91°S)0V39 
0-Х30М1 


51М312144302 VNOISN3WIG-NON QV3 


LIN20*HOIdON*9NVXVW* VL3H1V* L0Z* J410V* 3410. /XIS/ NOHHOO 
u1410*u*Isd*1130*8130*$130*0*v13HL /3A14/ NOWMWOO 

V3uvM vy* zv* Ly*gz*gA*gx* Lan* 10A * LGM* 102 

d' 1d0* 1O4* ZAI' AXI*ZXI* AI*XI* ZI* 92*9A* 9x*g* 14* 0H* q* W* 1* 34300N* 43301 
9W* 433023* 343022* 343024* 34300X* I2* I8* IV* IHd*0n*M* A*n. /4n04/. NOVIHOO 
1531 

103 “1НУО0531УЯМ“39У1158115711521151115915515015563: /ЗЗЯН1/ NOWWOD 
#21 

32“ 1“ Ххмимод“ 31v08*g9y1L* 9LL'SLL vLLIELLI2L* LL 29 * DI. /OML/ NOWWOD 
ХАН 

J1S* 1V8J1SSIVUOLLGL BL LA vL EL BN LN CONEA VA EN /INO/ NOWWOJ 
(0y)01 138 

(L11%009)M01S* (9*9)Y 8x 1438 

ZAI ZXI AXI ZI AI XI Sx WAU 

OLX 6A BA: LAIA SA VA EN ZALA Bx TWEE 
(9*9)yauviM* (€) r2* (6) 18* (€) 1v* (9)3* (9*9) vv* (22) 331 

3025N* (£2) 44302W* (£1) 333023* (22) 443002* (22) 343024 (91) 434300X 8x 1433 
(2-1)8x134* (H-4)8x* 1438 1191 1dHI 


“INdLNO JLIYM GNV S3ILI2OT3A 'SNOIIVN31390V “51М312133302 
-3'd ЗУЗМТ 31911979) “AZTIWNOISNAWIG-NON 9 3ZIlVILINI 
"SINIYA WILINT 8 SNOILIGNOD IVILINI 9 У1У0 0Уу334 --- 190034 NIVW 


O O Q 000 
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1130* L4vans * 313 vu*399un02*w9v11. (81*5)0v3y 
NOIIVWHOJNI 39NVH2 354002 0734 
LIN20* SuONO*HOIdON* J3JIGON (11*5)0y39 


НЭ114 9 Н1430 “358103 М04 58313НУНУд А11Л1115М35 0433 


1:191305313059813055134 (8159)0433 
ISd* IHd*V13HL (21*5)0v3y 

1qu* 100* 1dd*3*0*d (st*s)av3y 

10M* 10A* Lan*M*A*n. (st*s)avay 

L=1N9I 

(09 1=1*(1)01) — (666=0N3*p1*5)0V38 


'013 'S319NV *Q033dS 'Hid30 --- SNOILIQNOO TVILINI QV3H 


015563 (6155) 0У33 
8J LA IA GA (11 5)0V38 
21* p (£1*5)0v3y 


532\331$ 71041NO) dIHS 404 51МУ15М0Э 0934 


‘1M (EL S)OYJY 
ZAT*ZXI* AXI (215)av38 
ZI*AI*XI (21*$)av3y 
ОНЯ TEW An 
gz7*gA*ax (21's)av3y 
57*94*9x (21*s)av3y 


'013 'H19N31 ESSVW --- SITASIVILIVIVHO dIHS 0V3Y 


27417 (с1<6)0у38 
1*(1)19) (01*5)0v38 
1“(1)18) (01%6)0У38 
1*(1)1v) (01*5)av3y 


(еі 
(6-1 
(et 


чә чо < QOO 


с» С» с» 
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14), “Х9“, (sp) x45, (Q).tx^zts (aa) XL (Ld) ° XZ° ,(23S),°X€)1VW8M0J 2S 
(25*9)31 TM 
(,8130,*X4*,$130, X4* , 10M, * X8* , 10A, * X8* ,10nt 
6“ M XOLE A XOLE N, XOL Z XOLE A XOLE, X XL LI XS) 1VW03. L6 
( 16* 9)31IUM 
(/////*,8 
(930) ,*L't3*,- d3M011V 379NV 3AIG XVW,*x9*,(930) ,*t' 635,5 0340117 
V 319NV 3NV1ld Nu3ls ХУМ,“ Хог“ //“ ,(9308) .*t^vt4*,2 ОЗМОЛЛМ 319NV 3NV9 
14 М08 ХУМ, Х/15,(930) ,5t'vt4*,5 aanoTiv u3aaans xvW,*xoz*//*,(23s/S 
930) ,41781%,- ЗЈУМ 3NV1d NH318s,*xs*,(23s/930) ,°* LU 'b3*,= 31VY 3NV1b 
d MOg,°*Xt°,(523S/930) ,*t^t4*,2 31vu 30001, x02*//*,(930) «07945,6 
- 394002 0343090,*X01*,(13) ,*0'54*,2 1434 03353090, X€1* (23S) ,*Z 
L'£J*,- 1 vi130,*x02*//*,(93a0) 000 = 3$3109 9№113\1$,‘Х6‘ , (14) ,*01. 
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8662°{6/1Н4=1Н4 
6889" L«10M=10M 
688971х,10Л-10Л 
688971х100-100 
6889 L«MzM 
6889" |» Л=Л 
6889" | »П=П 
6889" | *0N=0N 
0300 
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0232382 . 


"00152582 
`O=IHdV 
`0=V13H1V 
`0=441@\/ 
`0=107 

: "0-10Х 
"0-10Х 
ши 
(00:0-Х 
`{=811 
15] 

" [591 
15-81 
“Ге? 
“БЕН 
12211 
“ЦІ 
"120 
“l=6L 
“1281 
[51] 
“1491 
[261 


OOO 








65 


100»113040-0 

і04х1130-аға 

19М» 1130+М=М 

° 10) 11304454 
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18 01 09 (140:0:0371)41 
86627/5/Ң2140М-Н2140М 
8562 "1G/SUINO=SYINO 
986562" LG/9NVXVW-ONVXVMH 
8562" LG /XWY31S=XWYILS 
8562 ° (5 /1\4315=1\4315 
8562" LG /XVWMOS-XVWMOG 
8962 ° 465/31Vu08-31V408 
86562" 1S/ Uv ani - Luvand 
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Nu31S TIVO ((HOLIdON*19'V13H1V) "QNV^ (LI34207 117 3410V) ) HI 
NY3LS 11У3 (1182071974410М)41 

налаза 7772 (331001' 19" 33107)31 

3399139 1772 (("0"314" 4738) "40" (542М10" 19" 2542) )41 
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(154-1582-14)598У0-2582 (:0:19715МЭ)41 
(154-1531002 )5890-2582 (:0:31:158Э)41 
Id-3Sun02-1592 

69191 °5=Та 

(IHd)Sava-IHdv 

(154)58У4-154У 

(VLIHL)SAVA=VL3HLY 

(4410) saVd=4JIQV 

Н14300-7-4410 
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10Ax11304A-A 

10X«11304X-X 
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(113a«(((ISd)NISQ«(V.13H2 
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input matrix of dimension N by N containing the coefficient 
matrix of the equation AX = B. | 

Оп output, A is replaced by the LU decomposition of a 
rowwise permutation of A. 

number of right-hand sides. (input) 

order of A and number of rows in B. (input) 

number of rows in the dimension statement for A and B in 
the calling program. (input) 

input matrix of dimension N by M containing right-hand 
sides of the equation AX = B. 

On output, the N by M solution X replaces B. 

input option. 

If IDGT is greater than O the elements of A and B are 
assumed to be correct to IDGT Decimal digits and the 
routine performs an accuracy test. 

If IDGT equals zero, the accuracy test is bypassed. 
work area of dimension greater than or equal to N. 
error parameter 

terminal error = 128 + N 

N = 1 indicates that A is algorithmically singular. 
Warning error = 32 +N. 

N = 2 indicates that the accuracy test failed. The 
computed solution may be in error by more than can be 


accounted for by tne uncertainty of the data. 


CALL LEQTIF (A, M, N, IA, B, IDGT, WKAREA, IER) 
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